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solute  attenuation  rates  due  to  water  vapor  and  moist 
nitrogen  have  been  measured  in  the  laboratory  at  138  GHz, 
282  and  300  K,  pressures  up  to  1.5  atm,  and  relative  humid¬ 
ities  from  80  to  100  percent.  The  computer-controlled 
measuring  system  is  comprised  of  a  millimeter  wave  reson¬ 
ance  spectrometer  (0.15  km  effective  path  length)  and  a 
humidity  simulator.  Several  shortcomings  of  earlier  meas¬ 
urement  attempts  have  been  rectified.  The  data  are  inter¬ 
preted  as  a  water  vapor  continuum  spectrum  consisting  of 
two  terms,  namely  strong  self-broadening  (H2O-H2O)  plus 
foreign-gas-broadening  (H2O-N2)  contributions.  Implica¬ 
tions  of  the  new  results  for  modeling  atmospheric  EHF 
window  transparencies  and  for  revising  established  H2O  line 
broadening  theory  are  discussed. 

Key  Words:  Atmospheric  radio  wave  propagation;  high- 

humidity  spectrometer;  millimeter  wavelength 
region;  water  vapor  continuum  absorption 
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Introduction 


Millimeter  wave  propagation  through  the  atmosphere  is 
affected  by  meteorological  phenomena  such  as  humidity,  fog, 
clouds,  and  precipitation.  System  applications  require  a 
detailed  knowledge  of  atmospheric  propagation  limitations 
and  their  dependency  upon  relevant  meteorological  parame¬ 
ters.  The  influence  of  the  intervening  atmosphere  upon 
radio  waves  can  be  quantified  by  a  propagation  model.  In 
clear  air,  water  vapor  and  oxygen  absorption  cause  fre¬ 
quency-dependent  signal  attenuation,  propagation  delay,  ray 
bending,  and  medium  noise.  Prediction  of  these  propagation 
effects  is  accomplished  by  a  broadband  model  of  the  complex 
refractivity  N  (1).  To  cover  frequencies  up  to  300  GHz,  a 
reliable  expression  for  N  requires  more  than  200  coeffi¬ 
cients.  The  refractivity  N  is  formulated  in  units  of  parts 
per  million  (ppm)  and  consists  of  three  components 

N  =  Nq  +  D(f)  +  jN"(f)  ppm 

where  NQ  is  the  frequency  independent  refractivity,  and 
D(f)  and  N"(f)  are  frequency-dependent  dispersion  and 
absorption  terms  arising  from  the  molecular  spectra  of  H20 
and  02. 

The  specific  rates  of  the  propagation  medium  are  power 
attenuation 


a  =  0.1820fN"(f)  dB/km  (1) 

and  delay 

8  =  3.336[Nq  +  D (f ) 3  ps/km  ,  (2) 

where  the  frequency  f  is  in  units  of  GHz. 

This  paper  addresses  millimeter-wave  water  vapor  at¬ 
tenuation  rates.  There  has  been  a  long-standing  discrep¬ 
ancy  between  theoretical  spectroscopic  predictions  and 
experimental  data  obtained  from  both  laboratory  and  field 
measurements  (3,4).  Marked  differences  occur  in  the  window 
ranges  located  around  35,  90,  140,  and  220  GHz  (2),  where 
experimental  attenuation  exceeds  the  theoretical  values. 
Since  excess  factors  between  2  and  5  are  involved,  it  has 
been  speculated  that  unidentified  absorbers  related  to  H20 


Atmospheric  Water  Vapor  Continuum  Mow  300GHz 


209 


(dimers,  clusters,  etc.)  might  be  involved.  Results  for  a 
from  controlled  laboratory  experiments  in  the  140  GHz 
window  are  reported  which  help  to  clarify  the  water  vapor 
excess  problem. 


Water  Vapor  Millimeter  Wave  Spectroscopy 

Water  in  both  vapor  and  liquid  states  obscures  the 
millimeter  wave  transparency  of  the  troposphere.  Frequency 
dependence  and  magnitude  of  HgO  attenuation  are  both 
distinctly  different  for  the  same  (1  mm)  absorber  thickness 
in  water,  rain,  suspended  droplet,  moist  air,  and  pure 
water  vapor  form,  as  illustrated  in  Figure  1  (2). 


Figure  1.  Normalized  specific  H2O  attenuation  over 
the  frequency  range  5  to  1000  GHz  for  five  phase  states: 
water  (ft  rain  (ft  suspended  hydrometeors  Q),  moist 
air  (ft  (numbers  shown  indicate  selected  H2O  attenua¬ 
tion  peaks) ,  and  pure  water  vapor 
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Water  vapor  attenuation  aw  depends  on  absolute  humidity 

v  =  7. 219e0  g/m3  (3) 

where  e  is  the  partial  vapor  pressure  in  kPa  =  10  mb  and 
0  =  300/T  is  a  relative  inverse  temperature  indicator. 

Water  vapor  pressure  e  is  limited  by  relative  humidity 

RH  =  41.51e0-5lO(9‘8340-1O)  £  100%  .  (4) 

The  attenuation  has  its  origin  in  the  absorption  lines 
of  water  vapor.  In  principle,  aw  can  be  evaluated  from  a 
line-by-line  summation  of  all  line  features,  each  contribut¬ 
ing  at  a  given  frequency  as  prescribed  by  a  line  shape 
function  F(f).  Six  reported  H2O  line  data  bases  (LB)  are 
indexed  in  Table  I. 

Table  I.  Available  data  bases  for  H2O  absorption  lines 


ID 

Frequency 

cutoff 

Number 
of  lines  (n) 

Remark 

Ref 

THz 

LBl 

537 

38350 

AFGL  Tape 

5 

LB  2 

126 

17201 

French  version 

6 

LB  3 

30 

2277 

Rotational  Spectr. 

7 

LB4 

3 

151 

Radio  Astronomy 

8 

LB5 

1 

30(+a  ) 

EHF  approximation 

1 

LB6 

1.1 

17  c 

Russian  version 

18 

The  extreme  low  frequency  wings  of  about  180  stronger  (i.e., 
at  vQ  more  than  105  times  the  value  of  a  at  100  GHz)  lines 
below  12  THz  make  contributions  in  the  EHF  range.  Six  dif¬ 
ferent  standard  line  shape  functions  F(f)  are  available  to 
assess  these  contributions.  Usually,  the  shape  functions 
are  governed  by  two  spectroscopic  parameters:  line  center 
frequency  v0  and  width  y.  Recently  a  third  parameter,  the 
population  "frequency", 

vB  =  2kT/hc  *  12510/0  GHz,  (5) 

has  been  introduced  (4b,  4h).  When  these  shape  functions 
are  normalized  to  unity  at  the  center  v  ,  the  functional 
dependences  given  in  Table  II  are  obtained. 
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Table  II.  Normalized  pressure-broadened  line  shape 

functions 


Normalization:  a  =  a  (v  )F(z)  for  a  single  line, 

woo 

a  =  Y/vo  ,  b  =  vo/vB  ,  z  =  f/vo 


ID*  F(z)  6xl03 

_  _  (Eg. 7) 


SL 

a2z  [a2+(l-z)2]_1 

1.9 

FL 

a2z  { [a2+(l-z)2]-1  -[a2+(l-z)2]-1} 

-0.001 

GR 

a2z2  fa2z2+(l-z2)2/4]-1  =  (FL) 

-0.001 

VW 

a2z2{[a2+a-z)2]-1  +[a2+(l+z)2)'1} 

637 

MGR 

(GR)z 

8.2 

MVW 

(VW) [tanh(bz)/tanh  b]/z 

16.6 

SL  =  short  Lorentzian,  FL  =  full  Lorentzian, 

GR  *  Gross,  VW  =  Van  Vleck-Weisskopf , 

MGR  =  modified  GR,  MVW  =  modified  VW  (10). 

The  intensity  distribution  of  a  strong  H2O  line  was 
traced  with  the  shape  functions  listed  in  Table  II  and 
the  example  is  exhibited  in  Figure  2.  Within  the  line  core, 
that  is  F(z)  10  or  z  -  0.8  to  1.2,  all  shapes  are  alike 
and  it  is  inconsequential  which  one  is  used.  On  the  other 
hand,  far-wing  intensities  depend  very  much  on  the  chosen 
frequency  distribution  function.  For  example,  the  relative 
contribution  of  the  line  shown  in  Figure  2  can  vary  at 
100  GHz  over  a  range  l(MGR) :75(SL) . 
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Figure  2.  Normalized  pressure-broadened  line  shapes 
F(z)  (see  Table  II)  applied  to  a  single  H2O  line. 

The  line  is  actually  centered  at  vD  =  1669.907  GHz 
with  a  peak  attenuation  rate  aQ  =  2.3  *  105  dB/km 
in  sea  level  air  (0  =  1,  v  =  10  g/m3). 

The  EHF  (30  to  300  GHz)  propagation  model  developed 
by  us  (1)  uses  for  reasons  of  computational  economy  30 
local  H2O  lines  below  1  THz  (LB5,  Table  I).  A  theoretical 
assessment  of  EHF  attenuation  a  =  was  accomplished 

on  the  basis  of  LB3  by  subtracting  the  line  contributions 
below  1  THz  from  all  2277.  The  far -wing  contributions 
from  pressure-broadened  lines  above  1  THz  are  then  fitted 
below  300  GHz  to  a  water  vapor  continuum  absorption 

«c  =  k(p/101)vX9y(f/100)u  dB/km,  (6) 

where  p  is  dry  air  pressure  in  kPa.  One  recognizes  ex¬ 
ponents  for  vapor  concentration  (x),  temperature  (y) ,  and 
frequency  (u)  dependences.  Applying  various  shape  func¬ 
tions  (Table  II)  to  LB3  led  to  the  results  summarized  in 
Table  III.  As  expected,  the  scale  factor  k  varies  over  a 
wide  range.  Also,  one  notices  that  all  three  exponents 
are  influenced  by  the  selected  shape  function. 
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Table  III.  Theoretical  EHF  water  vapor  continuum 
coefficients  (Eq.6) 


Line  Shape 

k 

x(v) 

y(T) 

u  (f  ) 

MGR 

0.0014 

1.00 

2.0 

2.5 

MVW.VW 

0.0064 

1.01 

1.4 

2.1 

GR.FL 

0.012 

1.02 

1.4 

2.0 

SL 

0.064 

1.05 

1.2 

1.2 

The  functions  were  tested  for  finite  integrated 
(f  =  0  to  “)  line  absorption  to  decide  which  shape  is  most 
plausible.  A  limited-range  numerical  integration  resulted 
in  residuals. 


L000 

=  (l/7ra)/F(z)dz  -  1, 


0/F< 
0.001 


(7) 


and  it  follows  (see  Table  II)  that  some  deviate  consider¬ 
ably  from  the  ideal  6=0.  Based  on  Eq.(7),  the  shapes  VW, 
MVW  and  MGR  can  not  be  considered  valid  intensity  descrip¬ 
tors  of  a  pressure-broadened  line  over  a  wide  frequency 
range  (22) — and  the  others  are  probably  overtaxed  by 
expecting  them  to  predict  eight  significant  orders  of 
attenuation  rate. 


Water  Vapor  Continuum  Absorption 


The  attenuation  rate  a  due  to  a  single  line  is  calcu¬ 
lated  with  Eq.(l),  where  the  product  of  strength  S  and 
shape  F  defines 

N"(f)  =  SF(f )  ppm  .  (8) 


The  line  properties  are  formulated  in  Ref.(l)  at  0  =  1 
(300  K)  by  a  set  of  spectroscopic  coefficients  bj  to  b3. 

In  the  low  frequency  far  wing  (y<<f<<v  ),  any  of  the  theo¬ 
retical  line  shapes  in  Table  II  is  proportional  to  the 
width  y;  however,  individual  continuum  shape  approximations 

F  (f  ->•  0)  = 
c 


4f2y/v£  [MGR],  2fy/v^  [MVW.VW], 


4fy/v^  [GR.FL],  y/v2  [  SL  ] , 

predict  successively  higher  contributions  (see  Fig. 
Table  III). 


(9) 

2, 
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Strength  S  and  width  y  of  a  typical  single  H2O  line,  pres¬ 
sure-broadened  by  a  mixture  of  water  vapor  (e)  and  dry  air 
(p) ,  are  given  by 


S  -  b^eO1*5 

kHz 

(10) 

y  -  b3(m*e  +  p)0 

GHz 

• 

(ID 

Water  vapor  continuum  attenuation  a  due  to  n  H2O  absorption 
lines  described  by  FL  or  GR  shapes  can  be  approximated  with 
Eqs.(l),  (8)  to  (11)  in  the  form 

a  *  0.73f2e(m*e  +  p)02*5  *  £  (b^/v3)  .  .  (12) 

c  i  =  1  °  1 

The  continuum  spectrum  consists  of  a  self -broadening  term 
(«  e2)  and  a  foreign-gas-broadening  term  («  ep) , 


a  =  (k  e  +  krep)0,y 
c  s  I 


dB/km 


(13) 


where  line  shape  theory  predicts  the  coefficients 

k  =  3. 5f2  l  =  mk 

s  i  =  1 

and  the  exponent 


y  =  2.5 


The  ratio 


(14) 


(15) 


m 


ks/kf 


=  4.80 


(16) 


is  the  pressure-broadening  efficiency  of  pure  water  vapor 
with  reference  to  dry  air,  determined  by  measurements  in 
line  cores  (15).  There  is  ample  evidence  (4a,  10,  12) 
that  in  the  EHF  range,  besides  local  lines,  the  low  fre¬ 
quency  wing  of  the  self-broadened  rotational  H2O  spectrum 
centered  around  9  THz  is  much  stronger  than  predicted  by  Eq. 
(12).  A  careful  review  of  H2O  continuum  absorption  by 
Burch  (10),  including  all  of  his  reported  data,  uncovered 
increases  of  the  broadening  efficiency  in  the  far  wings  of 
lines  ranging  between  m  -  20  at  300  GHz  and  m  -  50  at 
30  GHz  (see  later,  Table  IV). 


Atmospheric  Water  Vapor  Continuum  below  300GHz 


215 


The  laboratory  experiment  described  in  the  next  section 
was  prepared  to  measure  absolute  attenuation  rates  a  at 
three  EHF  window  frequencies  (i.e.,  92,  138,  and  220  GHz) 
under  conditions  reaching  water  vapor  saturation.  Ab¬ 
sorption  by  neighboring  lines  is  small,  so  that  to  first 
order  a  -  a  .  Attenuation  data  are  expected  to  behave  as 
predicted  by  Eq.(13).  The  ongoing  experimental  studies  are 
aimed  at  securing  accurate  coefficients  k  ,k^,  and  y,  as 
well  as  looking  for  possible  anomalies  not  supported  by  the 
pressure-broadened  rotational  spectrum  of  the  H2O  molecule. 
The  last  section  of  this  paper  presents  first  attenuation 
results  obtained  at  138  GHz  and  discusses  their  implica¬ 
tions  to  EHF  propagation  modeling. 


Millimeter  Wave  Spectrometer  for  High  Humidity  Studies 

Controlled  experiments  that  simulate  atmospheric 
conditions  provide  test  cases  for  studying  specific  contri¬ 
butions  to  N  in  isolation.  Assessments  of  basic  physical 
principles  underlying  the  attenuation  rate  a  are  difficult 
to  make  from  measurements  in  the  actual  atmosphere.  Re¬ 
liability,  precision,  and  scale  of  supporting  meteorologi¬ 
cal  data  compromise  the  quality  of  most  field  observations. 
The  expected  EHF  window  attenuation  is  low  in  a  spectro¬ 
scopic  sense  (<10  dB/km)  and  requires  long  path  lengths 
(>100  m)  for  a  reliable  detection.  A  semi-confocal  Fabry- 
Perot  resonator  with  a  high  Q  value  (437,000  at  138  GHz) 
offered  a  good  compromise  between  compactness  and  detection 
sensitivity.  The  computer-controlled  measuring  system 
consists  of  the  millimeter  wave  resonance  spectrometer  and 
a  humidity  simulator.  An  insulated  box  contains  a  high- 
vacuum  stainless  steel  vessel  which  houses  a  temperature- 
controlled  mini-lake  (10  cm  across)  and  the  radio  test  path 
(see  Figure  3).  Signals  at  frequencies  between  30  and  300 
GHz  can  be  employed  over  equivalent  free-space  path  lengths 
between  500  and  50  m;  temperatures  are  controlled  to  better 
than  1/100  of  a  degree  Celsius;  pressure  ranges  over  seven 
orders  of  magnitude  (102  to  10  5  kPa) ;  and  relative  humid¬ 
ity  is  varied  between  0  and  100  percent. 

Schematical  diagrams  of  physical  (Figure  3)  and  elec¬ 
tronic  (Figure  4)  arrangements  convey  an  overview  of  the 
experiment.  Quasi-static  diffusion  mixing  was  chosen  to 
generate  well  defined  humidities  inside  the  spectrometer 
cell.  A  temperature  controlled  water  reservoir  serves  as 
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s  Pressure  P  =  I0's  to  120  kPo 

*  Humidity  =  see  diagram  below; 


Figure  3.  Cross-section  of  millimeter 
wave  resonance  spectrometer  and  built-in 
humidity  generator.  Note  three  different, 
extremely  stable  temperature  zones,  Tj(0i) 
to  13(03).  Relative  humidity  inside  the 
cell  is  (Eq.  4):  log  (RH)  = 

2-9.834  (0]  -  ©2)  +  5  (logOj  -  log02)  %. 

Absolute  humidity  v  is  shown  in  the  diagram. 

the  vapor  source.  Electropolished  stainless  steel  was  used 
exclusively  as  construction  material.  Various  hydrophobic 
coatings  were  studied  as  possible  means  for  neutralizing 
the  absorption/ desorption  cycle  on  surfaces  exposed  to  water 
vapor  (4d),  but  were  abandoned  in  favor  of  slightly  heating 
the  mirrors.  Four  fast-responding  (t  <  Is)  temperature 
sensors  inside  the  cell  diagnose  any  disturbance  of  the  gas 
equilibrium.  Data  acquisition  is  computer-controlled. 
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Figure  4.  Schematic  of  high- humidity  millimeter  wave  spec¬ 
trometer  and  computer-controlled  data  acquisition  system. 


User-friendly  software  allows  the  choice  of  several  modes 
of  continuous  operation  over  time  scales  from  one  minute  to 
many  days.  The  resonance  signal  that  carries  the  absorp¬ 
tion  information  is  detected  with  a  digital  waveform  ana¬ 
lyzer.  The  analyzer  performs  the  following  tasks  under 
computer  control:  averaging  (1  to  100*),  baseline  normal¬ 
ization,  peak  value  (P-P)  and  bandwidth  (BW)  measurements. 


Detection  Principles 

The  dimensionless  complex  refractivity  N  is  a  macro¬ 
scopic  measure  of  the  interaction  between  gas  molecules  and 
mm-waves.  Real  and  imaginary  parts  are  related  by  the 
Kramers-Kronig  formula.  The  response  curve  of  a  gas-filled 
resonator  is 

a(f )  -  {Q_1  +  Qm_1  +  j2[(f/fR)(l  +  N’xio"6)-!]}-1,  (17) 
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where  Q  is  the  loaded  Q-value,  fR  is  the  resonance  frequency 
of  the  evacuated  resonator,  and  a  Q-value  of  the  gas  is  de¬ 
fined  by 

=  91021  x  f  /a  =  106/2N"  .  (18) 

M  K 

Two  modes  of  spectrometer  operation  follow  from  Eq.(17) 
when  the  evacuated  resonator  is  slowly  subjected  to  a  gas 
pressure  P; 


Absorption  mode.  The  resonator  is  swept  periodically 
across  its  resonance  fR  in  order  to  measure  a.  The 
shape  of  the  response  curve,  Eq.(17),  appears  as  a  pulse 
series  a(t)  when  rectified  with  the  detection  law  r.  The 
decrease  in  the  maximum  amplitude  a  (P  =  0)  at  fR  to 
values  a\ (?)  is  sensed  by  a  peak  defector.  The  attenuation 
rate  is  obtained  from 

a  =  4.343[(a  /ai)1/V  -1]/LD  dB/km  ,  (19) 

O  K 

where  the  effective  path  length  is  given  by 

Ld  =  4.771  x  10~5Q/fD  km  .  (20) 

K  K 

From  a(t) ,  the  peak  detector  records  an  absorption  pressure 
profile 

ax(P)  =  a  / { [a(P)Lp/4. 343]  +  1}B  ,  (21) 

0  tv 


where  the  total  pressure  P  =  e  +  p  is  measured  and  the 
expected  attenuation  rate  a  is  predicted  by  Eq.(13).  A 
digital  waveform  processor  (see  Figure  4)  resolves  amplitude 
changes  to  about  0.1  percent,  allowing  a  detection  sensi¬ 
tivity  of  (16) 


8  •  » 


a 

138  GHz, 


,  >  92fD/Q 

min  —  R 

Q  =  437000, 


a  . 
mm 


-  0.05  dB/km) . 


(22) 


Refraction  mode.  To  measure  N' ,  the  shift  in  resonance 
frequency  by  a  small  amount  to  f^  is  tracked  and  refraction 
is  given  by 

N'  =  106[ f R  -  fN(P)]/fR  PPm.  (23) 

Frequency  changes  of  fR  can  be  compensated  with  the  tuning 
micrometer,  allowing  a  resolution  of  =0.6  ppm. 


* 
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Specifications  for  the  138  GHz  Test  Series 

Resonator.  Fabry-Pdrot  ref lection- type,  semi-conf ocal 
arrangement,  10  cm  mirror  diameter,  20  cm  spacing;  micro¬ 
meter  tuning:  0.3175  mm/turn  with  1.3  x  10  6  mm  resolution; 
loaded  Q  value  at  f^  =  138.2  GHz:  Q  =  437,000;  effective 
path  length  (Eq.20):  L  =  0.151  km;  Fresnel  number:  6;  re¬ 
sponse  time  constant:  q/irf^  =  l.Olys;  coupling  coeffi¬ 
cient:  0.053;  coupling  hole:  circular,  0.65  mm  diameter, 
0.075  mm  double  Mylar  vacuum/ pres sure  seal. 

Excitation.  Linearily  frequency-modulated  klystron 
with  stabilized  center  frequency;  sweep  rate:  50  Hz,  sweep 
width  7  MHz  or  19600  us;  reflected  signal  off-resonance: 

150  mV. 

Detection.  Unbiased  Schottky  diode  (B  =  2)  with 
precision  attenuator  for  dB-calibration;  resonance  signal 
peak  ^sensitivity  (Eq.22):  ±  0.05  dB/km;  accuracy:  better 

than  ±  10  percent  for  a>l  dB/km. 

Data  Run.  An  example  of  measuring  water  vapor  attenua¬ 
tion  aw  by  detecting  the  resonance  peak  signal  ai (Pk)  with 
reference  to  vacuum  (a  )  is  shown  in  Figure  5.  Relative 
humidity  was  increased  from  RH  =  90  to  100  percent  over  a 
two-hour  period  by  raising  the  water  temperature  from  6.8 
to  8.4  C.  Condensing  water  vapor  added  several  tenths  of 
a  degree  to  the  cell  temperature;  however,  at  11  C  the 
mirrors  remained  free  of  condensation.  Some  of  the  data 
scatter  was  caused  by  drift  effects  of  the  spectrometer. 

The  calibration  of  a\  with  AT  (see  Figure  4)  determines  the 
dB-scale  and  is  completed  by  a  measurement  of  the  Q  value 
via  the  bandwidth  of  the  resonance. 


Results  and  Discussion 


Experimental  Results 

Pressure  scans  of  the  attenuation  rate  due  to  H2O 
absorption  were  made  at  the  frequency  fj  =  138.2  GHz. 

First,  pure  water  vapor  pressure  e  was  increased  to  a 
constant  relative  humidity  RH,  adjusted  for  values  between 
80  and  100  percent.  Then  nitrogen  was  injected  to  simulate 
moist  air,  although  broadens  water  vapor  lines  1.094 
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Figure  5.  Example  of  a  time  series  data  plot  from  five 
instruments,  reading  temperature:  Tj  =  mirror,  T^,4  = 
cell  center  and  bottom,  Qi  =  water;  and  peak  signal  Pk 
(water  vapor  absorption).  The  sequence  covers  a  two-hour 
period  during  which  water  vapor  saturation  was  achieved  by 
equalizing  the  temperature  of  8.4  C  between  cell  (T , 4) 
and  water  (Qi).  CAL  =  calibration  points  for  ai  (dB)  ,  BW 
(ms),  reference  power,  and  micrometer  setting. 


more  effectively  than  dry  air  (15).  The  test  frequency  fj 
falls  into  an  EHF  window  where  little  attenuation  by  local 
lines  is  expected  (see  Eq.27).  A  summary  of  reduced  and 
calibrated  data  is  shown  in  Figure  6.  A  data  fit  to 
Eq.(13)  at  300  K  amounted  to 

ct  =  0.197e^  +  0.00750ep  dB/km,  (24) 

resulting  in 

mX  =  26.3=  5.5m  .  (25) 

The  broadening  efficiency  m  is  determined  by  pressure 
pairs  ej  and  Pj — points  where  water  vapor  attenuation 
(a  e°)  and  moist  nitrogen  attenuation  [“  P(e-l)]  are  equal, 
i.e. 

mX  =  (Pi/ei)  -1 

Attenuation  at  points  A  (2.40  dB/km,  e  =  3.49  kPa)  and  B 
(5.10  dB/km,  P  =  106.7  kPa)  was  studied  with  special  care. 
Point  A  was  stable  over  a  40-hour  period  and  independent 
of  the  water  valve  position  (open/closed),  provided  a 
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Figure  6.  Pressure  dependence  of  specific  attenuation 
a  due  to  water  vapor  (e)  and  moist  nitrogen  (P)  mea¬ 
sured  at  138.2  GHz  and  relative  humidities  from  RH  =  80 
to  100  percent. 


slight  mirror  heating  (1W)  was  applied.  Increasing  the 
mirror  heating  to  5W  (+4°C)  did  not  change  the  result. 
Strong  anomalies  (condensation  effects),  on  the  other 
hand,  were  observed  when  the  mirror  heaters  were  off.  At 
point  B  a  very  long  mixing  time  (several  hours)  has  to  be 
allowed  when  the  water  valve  is  opened  against  110  kPa 
nitrogen.  Faster  mixing  occurs  when  nitrogen  is  blended 
with  water  vapor  already  filling  the  cell. 

Temperature  dependence  is  another  diagnostic  which 
can  be  used  to  support  the  theoretical  interpretation  of 
continuum  attenuation  in  terms  of  molecular  H2O  line 
broadening.  A  series  of  pure  water  vapor  pressure  scans 


222 


Urbr 


1 


i 


I 


✓ 


was  made  at  T  =  282  K  (see  Figure  5).  Temperature  limited 
the  saturation  pressure  to  1,102  kPa.  A  preliminary  fit 
of  these  data  yields  (compare  with  Eq.15) 

y  -  3.5  ±  1  (26) 

Pressure  dependence  indicated  by  Eq.(13)  has  been  ob¬ 
served  in  earlier  laboratory  (10  to  17,  23,  24)  and  field 
(4d,  4e)  experiments.  A  selected  review  of  laboratory 
data  is  given  in  Table  IV.  The  following  approximate  re¬ 
lations  can  be  deduced  from  this  data  body: 

-  kg(e2  +  ep/mx)  (in  between  lines); 
kg  *  f2  (excluding  the  67.7  GHz  data); 
mx  «  1/f^’^  (see  Ref.  10); 

y  -  5.5  (for  the  110  GHz  data). 


Table  IV.  Laboratory  results  for  water  vapor  attenuation 
a  fitted  to  Eq. (13)  yielding  the  coefficient  kg  and 
broadening  efficiency  m  . 


f 

kx 

s 

X 

m 

Broadening 

Agent 

T 

Ref . 

GHz 

dB/ (km 

kPa2) 

K 

31 

0.011 

49* 

Air 

318 

11 

67.7 

0.16 

44* 

Air 

283 

24 

110 

0.12 

35* 

Air 

304 

24 

110 

0.16 

35* 

Air 

292 

23 

110 

0.21 

35* 

Air 

274 

24 

138 

0.197 

26 

^3(29) 

(Air) 

300 

this 

138 

0.25 

- 

4 

282 

paper 

213 

0.62 

24* 

N2 

320 

4g»  17 

345 

2.4  ±  1, 

.5 

24 

N2 

300 

19 

411 

2.9  ±  0, 

.8 

18 

"2 

300 

19 

667 

4.1  ±  1. 

,1 

9.6 

*2 

300 

19 

883 

5.2  i  1, 

.4  9 

.4  ±  0. 

,9  N2 

300 

19 

Our  fit  to  reported  raw  data 
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Absolute  data  on  are  difficult  to  obtain.  In  spite  of 
existing  inconsistencies,  our  data  lend  support  to  the 
hypothesis  that  far-wings  of  self -broadened  H2O  lines 
centered  above  1  Tllz  arc  probably  the  major  source  for  the 
EHF  continuum  (10>.  Line  shape  theory  needs  to  be  revised 
to  account  for  the  minute  (<  10~5  the  value  at  line  center) 
contributions  in  the  extreme  low  frequency  wings  of 
pressure-broadened  lines  located  above  1  THz  in  the  rota¬ 
tional  H20  spectrum  (4b,  4h,  9,  21,  22). 

Water  Vapor  Continuum  and  EHF  Propagation  Modeling 

The  experimental  result  at  138.2  GHz  for  moist  nitro¬ 
gen  attenuation  (Eq.24)  has  been  incorporated  into  the  EHF 
propagation  model  that  is  based  upon  LB5  (1,25).  Reducing 
the  ep-term  by  9.4  percent  for  air-broadening  and  subtract¬ 
ing  line  contributions  due  to  LB5, 

ct£(LB5)  =  1.94  x  10-3  ep  +  0.0093  e2  dB/km,  (27) 

gives  the  pressure  dependence  of  continuum  attenuation  as 

c*c  =  =  4.86  x  10  3  ep  +  0.188e2  dB/km.  (28) 

Adding  frequency  and  temperature  dependences  to  Eq.28 
leads  to  a  moist  air  continuum  spectrum  for  LB5,  modeled  by 

N”  =  (CjepQ7  +  c2e20y*)f  x  10~6  ppm.  (29) 

The  coefficients  for  Eq.(29)  are  chosen  according  to  the 
following  scheme: 

c^  y  c2  y*  Remark 

Old  1.90  3.1  0  1  Gaut-Reifenstein 

New  1.40  2.5  54.1  3.5  Eqs.28,  15,  26 


The  "new"  continuum  spectrum  is  different  from  various  other 
empirical  forms  (summarized  in  Table  V)  which  have  been  pro¬ 
posed. 

Table  VI  compares  the  implications  of  the  old  and  the 
new  continuum  for  the  case  of  predicting  zenith  attenua¬ 
tion  A(dB)  through  a  model  atmosphere. 
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Table  V.  Experimental  water  vapor  continuum  coefficients 
(Eq.6)  reported  by  other  workers  (2). 


k 

x(v) 

y(T) 

u  (f  ) 

Local  Lines 

f,  GHz 

Ref. 

0.023 

1.00 

2.1 

2.50 

MVW 

50  to  500 

10 

0.039 

1.05 

3.0 

2.40 

GR 

20  to  350 

18 

0.049 

1.00 

2.1 

2.00 

VW 

15  to  500 

(*) 

0.061 

1.00 

2.1 

1.22 

GR 

150  to  410 

21 

Gaut-Reif enstein  model  (1). 


Table  VI.  Comparison  of  zenith  attenuation  A(dB)  predicted 
for  a  synthetic  atmosphere  (30  N,  July)  by  employing  the 
EHF  Propagation  Model  (1,25)  with  old  (Gaut-Reinfenstein) 
and  new  (Eq.29)  water  vapor  continuum  formulations. 


Ray  Path  Angle:  0.00  Degrees  from  Zenith 


Integrated  through  Model  Atmosphere:  pQ  =  101.35  kPa,  Tq  »  301.15  K, 
31  Height  Levels,  h  -  0  to  30  km. 


Integrated  Water  Vapor  V  *  52.69 


70.25  mm. 


Electric  Path  Length 


FREQUENCY (CHz) 

30.00 

35.00 

40.00 

45.00 

50.00 

55.00 

60.00 

65.00 

70.00 

75.00 

80.00 

85.00 

90.00 

95.00 

100.00 

105.00 

110.00 

115.00 

120.00 

125.00 

130.00 

135.00 

140.00 

145.00 

150.00 


«  2627.2 


ATTENUATION (dB ) 


.59 

.64 

.63 

.70 

.79 

.90 

1.13 

1.26 

2.13 

2.29 

28.16 

28.35 

148.98 

149.21 

23.87 

24.14 

3.19 

3.51 

2.42 

2.78 

2.36 

2.77 

2.48 

2.95 

2.69 

3.20 

2.94 

3.52 

3.23 

3.87 

3.58 

4.29 

4.07 

4.84 

5.47 

6.31 

11.89 

12.81 

5.38 

6.39 

5.51 

6.59 

5.92 

7.09 

6.43 

7.68 

7.03 

8.37 

7.75 

9.19 

old 

new 

2739.0  mm. 


ATTENUATION (dB) 


.72 

.88 

.74 

.96 

.91 

1.21 

1.27 

1.64 

2.29 

2.76 

28.36 

28.92 

149.09 

149.75 

24.16 

24.94 

3.49 

4.39 

2.75 

3.78 

2.72 

3.91 

2.89 

4.22 

3.14 

4.63 

3.44 

5.11 

3.79 

5.64 

4.20 

6.24 

4.75 

6.99 

6.23 

8.67 

12.73 

15.39 

6.31 

9.19 

6.52 

9.64 

7.03 

10.40 

7.66 

11.28 

8.41 

12.29 

9.32 

13.47 

old 

new 
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Summary 

Absolute  attenuation  rates  aw  have  been  measured  in 
the  laboratory  for  pure  water  vapor  and  for  mixtures  of 
vapor  and  nitrogen,  simulating  moist  air.  Initial  studies 
were  performed  at  138  GHz,  two  temperatures  (282  and 
300  K),  partial  vapor  pressures  up  to  saturation  (1.10  and 
3.48  kPa)  ,  and  total  pressures  up  to  160  kPa.  Preliminary 
results  were  acquired  for  pressure  and  temperature  depend¬ 
ences  of  o^.  The  data  behaved  in  a  fashion  similar  to 
that  expected  from  H20  line-broadening  theory  (Eq.13). 
Except  for  a  strong  self-broadening  component  (<*  e2), 
no  additional  anomalous  behavior  was  observed — even  at 
saturation  pressures  (RH  =  100%) — contrary  to  findings  by 
other  workers  (4,17,  etc.). 

In  our  case,  we  could  clearly  identify  anomalous  ab¬ 
sorption  behavior  (e.g. ,  unusually  high  rates,  extreme 
temperature  dependences,  hystereses  in  pressure  and  tem¬ 
perature  cycles)  as  being  instrumental,  that  is  condensa¬ 
tion  effects  on  the  millimeter  wave-active  parts  of  the 
spectrometer.  Considerable  effort  had  to  be  expended 
before  the  instrument  produced  consistent,  reproducible 
results. 

A  water  vapor  continuum  N"  was  formulated  for  the  EHF 
propagation  model  (1,25)  in  orSer  to  supplement  contribu¬ 
tions  from  30  H20  lines  below  1  THz.  At  138  GHz,  the 
local  lines  furnish  40  percent  to  air-broadening  but  only 
5  percent  to  self-broadening  of  the  water  vapor  attenua¬ 
tion  c^.  Practical  implications  of  the  continuum  to 
propagation  modeling  are  not  too  serious.  Up  to  humidi¬ 
ties  of  v  10  g/m3,  the  widely  used  empirical  Gaut- 
Reifenstein  continuum  (1)  is  roughly  confirmed.  At 
higher  values,  increases  over  the  GR  continuum  and  a 
nonlinear  dependence  on  absolute  humidity  are  predicted. 
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